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' We compute the ch spectrum from a first principle Salpeter equation obtained 
in a preceding paper. For comparison we report also the heavy-light quarko- 

^0 . nium spectrum and the hyperfine separations previously presented only in a 

\ graphical form. Notice that all results are parameter free. 
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^ , The spectrum and properties of the cb systems have been calculated various times in the 

Q-i! past in the framework of the heavy quarkonium theory [|I| . However the recent experimental 
{iLii observation of the B'^ meson P] has arisen new theoretical interest on the problem 

The mentioned spectrum has been considered again either from the potential and the lattice 
simulation point of view. A particular interesting quantity should be the hyperfine splitting 
that as for the cc case seams to be sensible to relativistic and subleading corrections in as. 
rS \ For the above reasons it seems to us worthwhile to present in this paper a calculation of 

\ the cb spectrum based on an effective mass operator with full relativistic kinematics which 
we have obtained in previous works and applied with a certain success to a fit of the entire 
quarkonium spectrum, heavy-heavy, light-light and light-heavy cases with the exception 
however of the cb case |§|J^. For comparison and completeness we report also numerical 
results for the light-heavy spectrum which we have given previously only in a graphical 
form. 

The mass operator was obtained by a three dimensional reduction of the qq Bethe- 
Salpeter equation introduced in P[. It has the quadratic form = Mq + [/, with a 



kinetic part Mq = wi + W2 = y mf + + y + and a "potential" that in terms of the 
istantaneous approximation of the B-S kernel is given by 



/iirni'\ ^ /wi + W2 f , Iw'i + w'2 , . 

kf^k =— — J- ^inst k,k J , ^ , (1) 

{271)'^ V 2WiW2 V 2WlW2 

k denoting the momentum of the quark in the centre of mass frame, and i = 1,2 the quark 
and the antiquark. 
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The B-S equation was derived from QCD first principles, taking advantage of the 
Feynman-Schwinger path integral representation for the "second order" quark propagator 
in an external field Q The only assumption used consisted in writing the logarithm of the 
Wilson loop correlator W = ^ {Ti P exp^^dx^^A^j)), as the sum of its perturbative expression 
and an area term 

i\nW = z(ln iy)pert + (rS^in, (2) 

a denoting the string tension. 

An explicit expression for U is given in Ref. 0. The perturbative part of such quantity 
was evaluated at the lowest order in ctg. However for as we have used the standard running 
expression 

-W) ^ (3) 

(with A^f = 4 and A = 200 MeV) cut at a maximum value as(0), to treat properly the infrared 
region 0. This amount to include important perturbative subleading contributions. 

Notice that, contrary to all the usual potential models, we have given the light quark 
current and not component masses in our treatment. Component masses of the usual order 
of magnitude can be recovered at a successive step as effective values in a semirelativistic 
reformulation Actually we have fixed such masses on typical values, = = 10 MeV, 
rus = 200 MeV, which are not adjusted in the fit (the results depend essentially on {k) 
and are very little affected by the precise value of the light quark masses). The other 
parameters of the theory are assumed as: = 1.394 GeV, mi, = 4.763 GeV, a = 0.2 GeV^, 
as{0) = 0.35. The first two are chosen in order to reproduce correctly the J/\l/ and the T(IS') 
masses, the string tension to give the correct slope for the Regge p trajectory, as(0) = 0.35 
to give the right J/'^ — rjc splitting. Notice that, consequently, the results reported in this 
paper are completely parameter free, with the exception of the cc{lS) hyperfine splitting. 

We have used in our calculations also the more conventional "linear mass" operator 
(or center of mass relativistic Hamiltonian) M = Mq + V (where V is defined hy U = 
MqV +V Mq + V"^) which makes easier a comparison with the usual phenomenological models. 
If we neglect the V"^ term, V is obtained from Eq. (|ip simply by the kinematical replacement 



V 2wiW2 V 2w{w'2 ~^ 4 ^WIW2W[W'2 ' 

This is the form we have used in Ref. @ (for some state however {V'^) is not negligible). In 
the calculations based on this linear formalism we have used the same values for the light 
quark masses as before, a fixed coupling constant as = 0.363 and taken rric = 1.40 GeV, 
rub = 4.81 GeV and a = 0.175 GeV^ 

Details on the numerical treatment of the eigenvalue equation are given in and [0. 



Second order propagator in the sense that it is defined by a second order differential equation; 
the quadratic form of the mass operator derives essentially from this fact. 
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In table | we have reported the cb spectrum as obtained by the quadratic and the hnear 
formahsm, together with the values presented in Refs. 0] and f^. The observed mass 
M(SJ = 6.40 ± 0.39 ± 0.13 GeV has to be referred to the 1 ^Sq state. For such state all 
calculations give very close results and reproduce equally well the experimental value within 
the errors. Larger discrepancies among the various methods occur for the excited states. 

In table || we have reported the spectrum for light-heavy mesons obtained by our for- 
malism in numerical form. We have considered the hyperfine structure but omitted the fine 
one. We have also reported the quantity Aavg defined as the average of the deviations of 
the theoretical values from the experimental data diminished by the experimental errors. 
Obviously Aavg provides a measure of the accuracy in reproducing the data and give an idea 
of the precision one can expect in the cb case. 

In table |ITI|, finally, we have reported the hyperfine splitting for the IS" and 25" states as 
obtained in the quadratic formalism and the Aavg quantity even for the channels for which 
we do not reproduce the results in full here. 

Notice the strong discrepancies with the data in the hyperfine splittings of the IS* light- 
light cases. This is obviously due to the chiral symmetry breaking problem and the related 
inadequacy of replacing the quark full propagator in the B-S equation with the free form, 
as implied in the three-dimensional reduction. For the rest, the agreement is good for the 
states involving light and c quarks, while the theoretical value tends to be too large for states 
involving b quarks. 

For comparison we can mention that in the linear formalism the hyperfine splitting turns 
out less good, being e.g. 97 MeV for cc{lS), 111 MeV for mc(15), 108 MeV for cs{lS). 
Likely such difference has to be ascribed to relativistic and subleading effects, taken into 
account in the quadratic formalism via Eq.@. 

In conclusion let us mention explicitly that Aavg as reported in table ^ do not include 
the states cc(4S') and bb{6S), (which are largely above threshold) and the 1 ^Sq and 1 ^Pi 
light-light states for the reasons recalled above. 
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TABLES 



TABLE L be quarkonium systems. Experimental mass equal to 6.40 ± 0.39 it 0.13 GeV. 



States 


quadratic formalism 


linear formalism 


Fulcher 


Lattice 




(GeV) 


(GeV) 


(GeV) 


(GeV) 


l^So 


6.258 


6.293 


6.286 


6.280 ± 0.200 




6.334 


6.355 


6.341 


6.321 ± 0.200 




6.841 


6.848 


6.882 


6.960 ± 0.080 




6.883 


6.881 


6.914 


6.990 ± 0.080 


S^So 


7.222 


7.221 






S^Si 


7.254 


7.245 






IP 


6.772 


6.762 


6.754 


6.764 ± 0.030 


2P 


7.154 


7.138 






ID 


7.043 


7.025 


7.028 




2D 


7.367 


7.346 
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TABLE II. Light-Heavy quarkonium systems. 


States 




experimental values linear formalism 


quadratic formalism 






(MeV) 


(MeV) 


(MeV) 


uc 












\ 


1869.3 ±0.5 1 
1864.5 ± 0.5 J 


1890 


1875 




' r>*(2oio)± 

r>* (2007)0 


2010.0 ± 0.5 1 
2006.7 ± 0.5 J 




2020 




D' 


2580 


2556 


2525 


2^Si 


D*' 


2637 ± 8 


2615 


2606 






2459 ± 4 ' 






IP 


£>|(2460)° 
* Di(2420)± 
L>i(2420)° 


2458.9 ± 2.0 

2427 ± 5 
2422.2 ±1.8 , 


> 2442 


2475 


^avg 






12 


21 


ub 










l^So 




5278.9 ± 1.8 1 
5279.2 ± 1.8 J 


OZoZ 


5273 


l^Si 


B* 


5324.8 ± 1.8 


5341 


5339 


2^So 






5878 


5893 


23Si 


B*' 


5906 ± 14 


5916 


5933 


IP 




5825 ± 14 


5777 


5792 


^avg 






34 


19 


SC 










I'So 


Dt 


1968.5 ± 0.6 


1999 


1982 


l=^Si 




2112.4 ± 0.7 


2107 


2120 


2^80 






2667 


2617 


2=^Si 






2729 


2698 


IP 


L>,j(2573)± 
D,i(2536)± 


2573.5 ± 1.7 
2535.35 ± 0.34 


i 9'i98 


2548 








21 


9 


sb 










1 ^So 




5369.3 ± 2.0 


5373 


5364 


l^Si 




5416.3 ± 3.3 


5433 


5429 


2IS0 






5974 


5985 


2^80 






6014 


6024 


IP 




5853 ± 15 


5848 


5859 


Aavg 






5 


4 
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TABLE III. Hyperfine splitting (MeV) 





UC 


sc 


ub 


oh 
SO 


cc 


CO 


hh 
00 


uu 


us 


ss 


IS 


145 


138 


66 


65 


115 


77 


86 


349 


298 


259 


Exp 


141(1) 


144(1) 


46(3) 


47(4) 


117(2) 






630.5(0.6) 


393.92(0.24) 


335.3(0.1) 


2S 


81 


81 


40 


39 


67 


42 


35 


135 


130 


127 


Exp 


57 








92(5) 






165(103) 






^avg 


21 


19 


9 


4 


20 




10 


19 


48 


18 
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